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Lack of inversion symmetry at a metallic surface can lead to an observable spin-orbit interaction. For certain 
metal surfaces, such as the Au(lll) surface, the experimentally observed spin-orbit coupling results in spin 
rotation lengths on the order of tens of nanometers, which is the typical length scale associated with quantum 
corral structures formed on metal surfaces. In this work, multiple scattering theory is used to calculate the local 
density of states (LDOS) of quantum corral structures comprised of nonmagnetic adatoms in the presence of 
spin-orbit coupling. Contrary to previous theoretical predictions, spin-orbit coupling induced modulations are 
observed in the theoretical LDOS, which should be observable using scanning tunneling microscopy. 

In the presence of time reversal symmetry [E{k,1) — E{—k,\,)] and spatial inversion symmetry [E{k,1) = E{—k,1)], 
no spin splitting can exist since E{k, f) = E{k, l). At a metal surface, however, spatial inversion symmetry is violated, and 
a spin splitting can therefore occur, i.e., E{k, t) ^ E{k, J,). The spin-orbit coupling in surface states was first observed by 
LaShell et al.' on the Au(lll) surface usin^hotoemission spectroscopy. The form of the spin-orbit interaction was found to 
be similar to the Rashba spin-orbit coupling^, which has been heavily studied in semiconductor heterostructures and quantum 
wells. Additional experimental-' and theoretical^ ^ ^ evidence have confirmed the presence of significant spin-orbit coupling 
on the Au(lll) surface. Although such a spin-splitting should, in principle, occur on all surfaces, the magnitude of the spin 
splitting depends very strongly on the nature of the surface. For instance, spin-orbit coupling has never been observed on either 
the Ag(ll 1) or the Cu(l 1 1) surfaces. This is due to the fact that the magnitude of the spin-orbit coupling is determined largely 
by the atomic spin-orbit coupling and the gradient of the surface state wave function at the nucleus^; theoretical calculations, 
which accurately predict the observed spin-orbit coupling on the Au(ll 1) surface, predict the spin-orbit coupling on the Ag(l 11) 
to be a factor of 20 smaller than the spin-orbit coupling on the Au(l 11) surfac^^El well outside the range of current experimental 
observation. In addition to the Au(lll) surface, photoemission experiments have discovered a variety of other metallic systems 
with spin-orbit coupling, such as on the Bi surfaces**, which exhibit an even larger spin-orbit coupling than that found on Au(l 1 1). 

Although most experimental observations of spin-orbit coupling in surface states are from photoemission spectroscopy, scan- 
ning tunneling microscopy (STM) has been used to observe spin-orbit interference in a magnetic sample^ and in nonmagnetic 
systems with very strong spin-orbit coupling, such as on the Bi(llO) surface'" and in Bi/Ag(lll) and Pb/Ag(lll) surface 
alloys". However, previous theoretical worlPhas argued that scanning tunneling microscopy (STM) could not be used to ob- 
serve the spin-orbit coupling in surface states; this argument was based on the assumption that the trajectories which interfere at 
the site of the STM tip are all one-dimensional in nature. Such trajectories do not undergo any net spin rotation, which results in 
the same standing wave pattern found in the absence of spin-orbit coupling. While the above argument is certainly true for the 
case of scattering from a single nonmagnetic adatom, trajectories involving multiple scatterers will undergo a net spin rotation, 
which will lead to spin-orbit induced modulations of the local density of states (LDOS), which should be observable using 

STM. .-.-^ 

Multiple scattering trajectories have been to be important in understanding the standing wave patterns observed 

in the LDOS for step edge^, for quantum corralJ^ formed by placin g adat oms atop a noble metal surface, and for quantum 
mirages generated by a magnetic adatom placed inside a quantum corraiSISl Previous experimental work has been conducted 
for quantum corrals on either the Cu(l 1 1) surface'^ or on the Ag(l 1 1) surfacJ^, where the neglect of spin-orbit coupling, as 
stated above, is completely justified^ ^ . However, this would not be the case for quantum corrals formed on the Au( 111) surface. 
In this work, multiple scattering theory in the presence of spin-orbit coupling is used to calculate the expected change in the 
LDOS for quantum corrals formed on surfaces with significant spin-orbit coupling, such as Au(lll). Numerical calculations 
performed for both a circular and a stadium quantum corral formed from nonmagnetic adatoms demonstrate that spin-orbit 
coupling can lead to observable changes in the LDOS. Understanding the effects of spin-orbit induced interference on metal 
surfaces will be important if such systems are to be used for future spintronics applications. 

The effective Hamiltonian for a surface state in the presence of the Rashba spin-orbit interaction is given by: 

where m* is the effective mass, a is the spin-orbit coupling strength, and Eq is an energy offset arising from the confine- 
ment of the electron to the surface. The eigenstates of H with energy E are given by ^'i(r) = exp(iA:i(i?) -7^1 + {4>)) 
and ^'2(7^ = exp(ifc2(-£') • r)\ — (0)), where the spin quantization axis for | ± {(j))) depends upon the momentum vector, 

k^2)iE) = k^2)iE) (cos(0)x + sin(0)y), where ki{E) = fcso + HE) and k2{E) = -k^o + k{E), with fcso = ^ and 
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k{E) = y (fcso)^ + ^™ ^^-^a) (foj- convenience, the energy dependence of fci, k2, and fc will not be explicitly written from 

now on). For a given value of 0, the spin states are | ± (</))) — (\/2)^^ (l+)2 i cxp(— i</))|— where li)^ are eigenstates of 

CT^. Due to spin-orbit coupling, the dispersion relation, E{k) = ^2m* "F '^1^1 + consists of two parabolic bands centered 

about ibfcso with the bottom of the bands occurring at an energy Eq — Esq (where Esq = .^Jf ) instead of at energy Eq. The 
dispersion relation is plotted in Figure [T|A) for kx — 0, where the spin states are | + (0)) = \+)x for the band centered at 
ky = fcso and | — (0)) = \—)x for the band centered at ky — — fcso- The full two-dimensional dispersion curve in the kx — ky 

plane can be found by simply rotating the dispersion curve in Fig. jlj^) using cxp (^—iO/hLz^ cxp {—i6/2az), where 9 e [0, tt) 

and Lz is the z-component of the angular momentum operator. 

In an STM experiment''', the bias voltage between the tip and the surface, V, can be changed in order to probe the local density 
of states at an energy Ep + eV (where Ep is the Fermi energy of the metal) by measuring the local conductance, ^ {Ep + eV, f), 
since l(,{Ep + eV,f) cx LDOS{Ep + eV,r) where LDOSiE,f) = 'Eg^±J2^{r',q\'i'^){'i'i,\r,q)S{E - E^). Thus in order 
to calculate the STM image, the LDOS{Ep + eV, r) must be determined. One method of determining the LDOS{Ep + eV, r) 
is by calculating the Green's function, G±{f 1,7^2, Ep + eV), and using the following relationship: 

LDOS{Ef + eV,r) -^Trace,spin [G+(f, f, Ep + eV) - G-{f, f, Ep + eV)] (2) 

Thus knowledge of the Green's function can be used to calculate the expected STM signal. 

The free-particle Green's function in the presence of Rashba spin-orbit coupling and for E > Eq — Eso is given bjl^SED 

r" r F) - Ti V G±iE,ri2) ±iexp{ie^2)Gg{E,n2)\ 

where = \fi - f2|, exp(±i6'i2) = [{fi - r2) ■ y ± i(^i - ^2) • x]/ri2, G^{E,ri2) = kiH^{kiri2) + k2H^{k2ri2), 
and G%{E,ri2) = kiH^{kiri2) - k2H^{k2ri2). Note that fcz < for energies Eo — E^o < Ep + eV < Eq, so that in 
this energy range, i?^(fc2?'i2) = {—^)"^^ {\k2\r12) in Eq. ([sjl. This results in a change in the LDOS{Ep + eV,r) when 
Eq — Eso < Ep + eV < Eq: for Ep + eV > Eq, the free particle LDOSftee{Ep + eV, r) is independent of energy and is 
given by LDOSfree{Ep + eV,f) = m* /{ttTl^), whereas for Eq - Eso < Ep + eV < Eq, the LDOS{Ep + eV, r) is dependent 
upon Ep + eV and is given by LDOSfieeiEp + eV, f) = m* /(jrh^) x kso/k. This change in the LDOSf^ee has been recently 
reported for STM measurements on Bi/Ag(l 1 1) and Pb/Ag(l 11) surface alloys" . 

In order to gain more physical insight into the transport between ri and r2, Gj.{fi,'r2,E) can be rewritten in terms of a 
complex amplitude multiplied by a "complex" rotation: 

Gl{n,r2,E) = Ti^bb^{E,n2)R(^^,±z±{E,n2),-^-f^ (4) 

where bb±{E,ri2) = ^J[G^{E,r 12)]'^ + [G^{E,ri2)]^, R{a,(3,"f) = exp(iaCTz) exp(i/3CTx) exp(i7CT2) is an ar- 
bitrary rotation operator with Euler angles a, (3, and 7, and z±{E,ri2) is a complex angle which is defined by 
cos[2;±(£;,ri2)] = G^{E,ri2)/bb±{E,ri2) and sin[z±(£', ria)] = G%{E, r 12) /bb±{E,r 12). Note that for a tra- 
jectory which goes from r2 to ri and then back to r2, no net spin rotation occurs, since G'j.(r2, n.)Gj.(ri, r2) oc 
{bb±{E,r^2)fR{'f,±z±{E,ri2),-'f)R{'f.Tz±{E,r^2),-'f) = {bb±{E,r^2)fl. 

In the presence of multiple adatoms, the total Green's function can be significantly altered from G^{ri,r2,E) due to the 
interference between the various multiple scattering trajectories. The total Green's function in the presence of N nonmagnetic 
adatoms/scatterers can be approximated as: 

JV 

G±{fur2.E) ^ Gi(rl,r2,S) + 5]Gi(rl,f,,S)?±G±(r;-,r2,-B) (5) 

where is the "s"-wave scattering amplitude, which is given by = h'^ lm*[eyip{±2\5j{E)) - 1]1, with 5j{E) being the 
scattering phase shift. In writing Eq. (jsjl, the scattering length of each adatom was assumed to be much smaller than 2tt /k 
(justifying the "s"-wave approximation) and the spin rotation length, n/kso, which allows one to associate the same scattering 
amplitude for both the ki and ^2 scattered waves (see for example Eqs. (32)-(33) of Ref.lSl). The unknown values of the Green's 
function at each scatterer n, G±(r„, r2, E), can be found by setting ri = r„ to give: 



G±(r;,r2,^) = Gl{f„,f2,E) + Y,Glirn,rj,E)sfG±{f,,r2,E) (6) 
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This resuhs in a system of equations which can be solved via a simple matrix inversion. With knowledge of G±(r„, r2, E) 
for each scatterer n, the total Green's function, G±(fi, r2, E) in Eq. ob, is determined, thus determining the LDOS by using 
Eq. ^. 

Consider first the simple case of a single nonmagnetic adatom placed atop a metal surface at Tj. The total Green's function in 
this case is given by: 

G± (fi , f2 , E) = Gl in , r2 , E) + sf G° {n , Tj , E) G° (r, , , E) (7) 
which results in a change in the LDOS{E,r) of ALDOS{E, t^q) = LDOS{E,ra) - LDOSfr^^{E,ra) = 



sin(2fcrnj)Im[s+]) (8) 



sin(2A;roj)Re[s+]) (9) 

for Eq — Eso < Ep + eV < Eq. Therefore, there exists a change in the ALDOS when Eq — &p < Ep + eV < Eq due 
to spin-orbit coupling, which is similar to the change observed in the LDOSf^e described earlieiE]. For the case of a single 
nonmagnetic adatom, this change in the ALDOS would be the only way to detect the presence of spin-orbit coupling, since 
the period of the spatial modulation in the ALDOS, 2k, can only be used to determine the effective energy of the surface state 

electron, k = ^ £cff "^ Spin-orbit coupling only shifts the effective bottom of the band from Eq to Eq — E^q, so measurement 

of k cannot, by itself, help to determine the presence or absence of spin-orbit coupling. The physical reason why spin-orbit 
coupling doesn't affect the LDOS in the presence of a single adatom is that for single scattering paths returning to the STM 
tip, no net spin rotation can occur, as shown in Fig. [TjB). This was the reasoning used to argue that STM couldn't be used to 
observe spin-orbit coupling for a surface state^. However, in the presence of multiple adatoms, mutliple scattering trajectories 
can generate a net spin rotation (Fig. [TjB))^ which allows the spin-orbit coupling to affect the LDOS in a nontrivial manner As 
mentioned earlier, such multiple scatteri ng trajectories are important in understanding the observed LDOS in quantum corrals 
formed atop noble metal surfacei^^EEll 

For the calculation of the L_D05 on the Au(l 1 1) surface, the following parameters were used': m* — 0.26TOe and — i?o = 
0.41 eV (Refl^^l), a spin-orbit coupling constant of a = 4 x lO^'^'^ eV-m (which is 10% smaller than the value given in Ref!^ 
and 21% larger than the value given in Ref.^l). These parameters give a Fermi wavelength of Af — 2it /k — 37. 4A and a spin 
rotation length of n/kso ~ 230. SA. It should be noted that this spin rotation length is about an order of magnitude smaller than 
the attainable spin-rotation lengths in semiconductor heterostructures, which is mainly attributable to the larger effective mass 
of the surface state electrons. 

In the following calculations, all adatoms were modeled as "black-dots'''^ where S{E) = ioo due to inelastic scat- 
tering of electrons into the bulk^^ (modifications of the theory for treating the adatom scattering as purely elastic^'' can 
also be performed too). In the simulations, each nonmagnetic adatom was placed on a hexagonal lattice at a position 
r= I ((6i + 62) 2? + {bi — i>2)y), where a = 2. 885 A for Au(ll 1), and bi and 62 are integers chosen to minimize |r — rd|, 
where is the desired location for each adatom. It should be mentioned that a hexagonal lattice is a simplified model of the 
actual Au(lll) surface, which undergoes a herringbone reconstruction^^. Such a reconstruction acts like a superlattice for the 
surface state electrons and modifies the electron density; however, such a reconstruction has minimal effect on the spin-orbit 
coupling as has been demonstrated by theoretical calculation J^l^ and is not considered in the following simulations. 

In order to illustrate the effect of spin-orbit coupling on the resulting ALDOS{E,r), simulations with and without spin- 
orbit coupling were performed at slightly different applied voltages but with the same effective energy, E^ff, in order that both 
simulations gave the same period in the spatial oscillation of the ALDOS{Ef.ff, f) in the presence of a single adatom, 2k = 

2\J 2m* Egff / . For example, if V was the applied voltage used in the simulation in the absence of spin-orbit coupling, then 
the applied voltage in the presence of spin-orbit coupling, V', would be given by eV' — eV — E50 = eV — 2.7 meV, with 
i?eff Ep — Eq + eV ~ Ep — Eq + eV + Esq. In order to consider only the contributions of spin-orbit to the ALDOS arising 
from multiple scattering trajectories, effective energies, E'eff > Eso, were only considered in order to avoid the intrinsic change 
in the ALDOS when < E^ff < Esq. Note that for the case of the Au(lll) surface, this intrinsic change in the ALDOS 
should in any case be unobservable since E50 = 2.7 meV is much smaller than the lifetime broadening" of 18 meV. 

Simulations were first performed on a circular quantum corral of radius 88. 7 A comprised of sixty nonmagnetic adatoms 
placed atop a hexagonal lattice The calculated ALDOS{Eef{, 0) at the center of the corral is shown in Figure|2]^A) as a function 



^ (^2") [-^t^b+iE, roj))^] , which, for kvQj ^ 1 can be approximated as: 

ALDOSiEF + eV,ro) = ( ^] (cos(2fcro,)Re[s^ 

k rQj VT?! J 

for Ep + eV > Eq, and as: 



2 

m 



ALDOS{Ep + eV,fQ) = I ^ j (cos(2fcroj)Im[s+ 

k Toj \ 
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of _Eeff in the presence (solid curve) and in the absence (dashed curve) of spin-orbit coupling. The ALDOS{Efif{, 0) without 
spin-orbit coupling has been shifted down for convenience. A very simple "particle in a box" modeP^can be used to interpret 
the ALDOS{E^s, 0) in Fig. ^A): in the absence of spin-orbit coupling and treating the quantum corral as a circular billiard 
with radius R ~ 88. 7A, the peaks in the ALDOS{Eeff, 0) mainly occur when i?eff is equal to an eigenenergy of the circular 

billiard, i?eff = En — ^ 2m*°'' where fc„ o is given by the solution to Jo{knfiR) = 0. This simple model predicts the peak 
locations in the ALDOS{E^^, 0) to within 10 meV for the first four peaks shown in Fig. |2|A). 

A similar model can be applied to the case of a circular billiard with spin-orbit coupling. In this case the eigenstates can be 
written as: 



i{r) oc exp(im0) 



which have an effective energy (shifted by —Esq for comparison to the simulations without spin-orbit coupling) given by 
Eeff = ?i^(fc"'™ + ^2 '™)^/ (8to*), which is determined by the condition: 

j^^iCfc^^'^i?) J„(fc2"'"i?) + J,n{k'r'R)Jm~i{k^''^R) = (10) 

The solutions to Eq. ( [To| i which can have nonzero amplitude at the center of the circular billiard, the degenerate states ^„ o(?^ 
and i(r), essential come in two types of eigenstates. The first type occurs at energies i?^ which are only about one to two 
meV smaller in energy than for the eigenstates Jo(fcn,oH) in the absence of spin-orbit coupling. These states, although pos- 
sessing some Ji(fc|7^) exp(±i6') character, are mostly Jo(A:|r|) in character, which leads to large peaks in the ALDOS{Ess, 0) 
at slightly lower energies than the corresponding peaks in the absence of spin-orbit coupling. The second type of eigenstate 



determined by Eq. ( 10 1 occurs at energies in between the aforementioned energies. These eigenstates, which are closely related 
to the m = ±1 eigenstates in the absence of spin-orbit coupling, J±i(fc„|7^) cxp(±i(?), possess a small amount of Jo(A;|r|) 
character due to spin-orbit coupling, which can lead to new, albeit small, peaks in the ALDOS{Ef.f(, 0) at these energies. The 
small peak in the ALDOS{E^f[ = 31meV, 0) (and more clearly shown in the inset in Fig. [2]^A)) corresponds roughly to such 
an eigenstate, which, for the circular billiard with spin-orbit coupling, has an energy of E^ff = 29.1 meV. 

Besides the small shift in the peaks of the ALDOS{Egf(, 0) and the small peak at E^g = 31 meV, the observed difference 
in the ALDOS{Esfi^ 0) with and without spin-orbit coupling is relatively small. However, the ALDOS at other places inside 
the quantum corral can show considerable differences when spin-orbit coupling is included. A slice of the ALDOS{Egff, r) 
through the quantum corral is shown in Figs. and|2|C) for the (C) second peak in ALDOS^Ees, 0) [i?eff — 58.7 meV 
(solid curve) and ii^eff = 60.4 meV (dashed curve)]and for the (B) third peak in ALDOS{Eef{,0) [iSeff = 144.1 meV (soHd 
curve) and E^ff = 145.8 meV (dashed curve)], where the black rectangles centered at ±88. TA correspond to the positions of 
the adatoms in the slice. Note that in the simulations, the ALDOS{Ei.g, f) is never calculated within 6A of the adatoms. As 
the electron bounces around in the corral, it undergoes an effective spin rotation due to spin-orbit coupling, which modulates 
the interference patterns seen in the quantum corral, resulting in an enhancement (Fig. |2|C)) or a decrease (Fig. ^B)) in the 
ALDOS{E^f(, r) near the edges of the corral. 

Besides the circular corral, another corral simulated in this work was a 78 adatom stadium billiard of dimensions 14lA by 
285 A, where the adatoms were again placed atop a hexagonal lattice. Figure [3] gives the ALDOS {Eef[,r) with and without 
spin-orbit coupling for roughly zero bias voltage between the tip and surface, i.e., i?eff = 410 meV. Calculations performed 
at different E^ff gave similar results (data not shown). As for the circular corrals, the ALDOS{Egff, r) was artificially set to 
zero within 6A of each adatom, which makes the adatom positions clearly visible in Fig. [sj^A). Although the general structure 
of the ALDOS{Esff, f) with and without spin-orbit coupling appears similar, spin-orbit coupling causes additional structure, 
such as splittings and intensity variations, to appear in the ALDOS{Eg{f, r) as is shown in Fig. |3jA). Such spin-orbit induced 
interference effects can be more clearly seen in Fig. |3jB), which plots a slice of the ALDOS{Egff, r) through the center 
of the stadium corral along the long dimension of the corral. As with the circular corral, changes in the amplitude of the 
ALDOS{Eeff, r) are seen near the adatoms (black rectangles in Fig. [sj^B)). However, spin-orbit coupling causes a splitting in 
the ALDOS{Egff, r) near the center of the stadium, where the peak to peak distance is roughly 18A. Additional splittings and 
modulations of ALDOS{Egff, r) can also be seen in Fig. [3|A). These calculations clearly demonstrate that spin-orbit coupling 
can generate significant changes to the ALDOS{E^ff^ r) in quantum corrals. 

In this work, we have examined the effects of spin-orbit coupling on the local density of states for quantum corrals formed 
atop the Au(lll) surface. Changes in the LDOS{Esf{, r) in both circular and stadium corrals indicate that spin-orbit induced 
interference effects should be visible using STM on the Au(lll) surface, contrary to previous theoretical arguments^. The 
modulations in the LDOS{Ef.f[, f) were a result of non-collinear multiple scattering trajectories, such as those found in quantum 
corrals, which can generate an effective spin-rotation in the presence of spin-orbit coupling. Sin ce th e previous experimental 
data on quantum corrals is quite good and can be accurately described by multiple scattering theorjl^^lil) jjjg predicted spin-orbit 
induced interference in such systems should also be experimentally observable. Furthermore, this work can also be extended 
to the case of quantum corrals comprised of magnetic adatoms, where, through the spin-orbit coupling of the surface state 
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FIG. 1: (A) The dispersion curve projected along fex = in the presence of spin-orbit coupling. The normal parabolic dispersion relation has 
been split into two parabolic curves, centered about ky ~ ±fcso, with the band edge occurring at an energy of Eq — Eso- Note that the spin 
state is | + (0)) = \+)x for the parabolic band centered at ky = +fcso and j — (0)) = \—)x for the parabolic band centered at ky = — fcso. 
(B) For a single scattering trajectory, spin-orbit coupling cannot generate a net spin rotation since R{9, /3, —0)R{9, — /3, —6) = 1. However, 
for non-coUinear multiple scattering trajectories, a net spin rotation can occur since R{6, ±/9, —9)R{5, ±7, —5)R{C^, ±A, —() 7^ 1. 



electronJ^, effective interactions between the magnetic adatoms can be generated. The methodology in this paper can also be 
used to calculate the LDOS for superlattices formed from localized structures. In the future, ab initio calculations of STM 
images in quantum corrals^^ with spin-orbit coupling will be performed, in addition to examining the effect of the herringbone 
reconstruction on the resulting LDOS in quantum corrals on Au(l 1 1). 
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FIG. 2: The ALDOS{E, r) for a circular quantum of radius 88. 7A comprised of 60 nonmagnetic adatoms, wliich are modeled as "black-dot" 
scatterers (S{E) = ioo). In (A), the ALDOS^E^n, 0) is plotted at the center of the coiTal as a function of the effective energy, E^g, with (solid 
curve) and without (dashed curve) spin-orbit coupling. The dashed curve has been shifted downward for convenience. Note that E^f can be 
converted into a bias voltage by using either —eF = i5F — i?o — -E'cff (dashed curve) or— el/' = Ef — Eo + Esq — E,:{i = Ef — Eq — Ei:ff + 2.7 
meV (solid curve), where i5 J? — iSo = 410 meV for the Au(1 1 1) surface. The following parameters were used in the simulation: m* — 0.26me 
and a = 4 X 10"^^ eV-m. The peaks in the ALDOS{E^fi, 0) occur when E^ft roughly corresponds to eigenergy for a circular billiard with 
(solid curve) and without (dashed curve) spin-orbit coupling. A small peak (shown in the inset) in the ALDOS{Ecff, 0) at Eeff ~ 31 meV, 
corresponds to an eigenstate of the circular billiard with spin-orbit coupling which is mostly Ji (fc|r|) exp(±i6') in character, but, due to 
spin-orbit coupling, does contain some Jo{k\r\) character, which can contribute to the ALDOS'(i5eff, 0). In (B) and (C), profiles of the 
ALDOS{Ec[[, r) through the quantum corral (the adatoms are denoted by the black rectangles) with (solid curve) and without (dashed curve) 
spin-orbit coupling for (C) the second main peak in the ALDOS{Es,ff, 0) [iJeff = 60.4 meV (without spin-orbit coupling) and i?eff = 58.7 
meV (with spin-orbit coupling)] and for (B) the third main peak in the ALDOS{Ecff, 0) [i?cff = 145.8 meV (without spin-orbit coupling) and 
Scff = 144.1 meV (with spin-orbit coupling)]. In both cases, substantial differences in the intensity of the ALDOS{Eaf, r) are observed, 
where the presence of spin-orbit coupling can either enhance the ALDOS{Ecf[, r) (Fig.|2|C)) or decrease the ALDOS{Es{f, r) (Fig. |2JB)). 
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FIG. 3: (Color online) Simulation of the ALDOS{Es[f, r) for a 78 adatom quantum corral stadium billiard of width 14lA and length 285A at 
Ecff = 410 meV, with (left) and without (right) spin-orbit coupling. Although the general features are similar, inclusion of spin-orbit coupling 
can enhance or diminish features in the ALDOS along with introducing additional peaks in the ALDOS. This can be more clearly seen in 
(B), where a slice through the center of the stadium along its long dimension has been plotted with spin-orbit (solid curve) and without (purple 
dashed curve) spin-orbit coupling. The black rectangles indicate the locations of the adatoms through the slice. Besides differences in peak 
intensity, a splitting of the ALDOS occurs at the center of the stadium with spin-orbit coupling (peak to peak distance of ~ ISA, which is 
absent when spin-orbit coupling is not included. 



